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in: 



We are reporting on the ongoing effort of the Monte Carlo (MC) modelling of NLO DGLAP QCD evolution in 
the fully unintegrated form. The resulting parton shower MC is performing on its own the NLO QCD evolution, 
contrary to all known programs of this kind which are limited to LO level only. We overview this new MC scheme, 
for the non-singlet subset of the gluonstrahlung diagrams. Precision numerical test of this new scheme are also 
shown. 
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1. Introduction 

In the so called factorization theorems, see for 
instance |1I2I3] , the scattering process with a sin- 
gle large transverse momentum scale (short dis- 
tance interaction), can be described in the per- 
turbative QCD |4|5|6j (pQCD) in terms of the 
on-shell hard process matrix element convoluted 
with the ladder part. The hard process is cal- 
culable up to a fixed perturbative order. The 
ladder part, calculated for each coloured par- 
ton entering (exiting) the hard process, is conve- 
niently described as a tree-like stream of partons, 
parametrized by the inclusive parton distribution 
function, PDF. 

The logarithmic response of the inclusive PDF 
to the change of the large kx scale is the so called 
DGLAP [7] evolution of the PDFs. This evolution 
was mastered for the inclusive PDFs up to NLO 
level in the early 80's, see for instance [8 9 , and 
recently even to the NNLO level [TU] . 

Alternatively the parton tree encapsulated in 
the inclusive PDF can be modelled using direct 
stochastic simulation of their four-momenta and 
other attributes, the so-called Monte Carlo (MC) 
technique. Here, the breakthrough was made in 
mid-80's, see refs. |11I12) . where the LO ladder 
was implemented in the so-called parton shower 
(PS) MCs. 



Standard LO level PSMC implements also the 
hadronization of quark and gluon partons into 
hadrons and is a workhorse in the software in 
all collider experiments. The advent of the LHC 
puts a challenging requirement on the quality of 
the pQCD calculations needed for the experimen- 
tal data analysis. Some of them will soon require 
that the PSMC is upgraded to the complete NLO 
level. This, however is not an easy task, mainly 
because classic factorization theorems [11213] are 
not well suited for the exclusive MC implementa- 
tion, but rather for inclusive PDFs. 

In this contribution we would like to report on 
the serious and successful attempt of solving the 
above problem of constructing NLO PSMC. We 
shall overview main technical points and present 
some numerical tests of the scheme. 

Possible profits from NLO PSMC include: 
(a) Complete set of "unitegrated soft countert- 
erms" for combining hard process ME at NNLO 
with NLO PSMC; (b) Natural extensions towards 
BFKL/CCFM at low x; (c) Better modelling 
of low scale phenomena, Q < lQGeV, quark 
thresholds, primordial k T , underlying event, etc.; 
(d) Porting information on the parton distribu- 
tions from DIS (HERA) to W/Z/DY (LHC) in 
the MC itself, instead in the collinear PDFs (uni- 
versality must be preserved) and more. 

The MC modelling of NLO DGLAP is not so 



much the aim itself - it will be rather a start- 
ing platform for many interesting developments 
in pQCD in many directions. 

Presently, we concentrate on constructing NLO 
PSMC for QCD initial state radiation (ISR) just 
for one initial parton such that: (a) it is based on 
the collinear factorization theorems |l|2|3j as rig- 
orously as possible, (b) we take scheme of Curci- 
Furmanski-Petronzio (CFP) as a main reference 
and guide (axial gauge, MS dimensional regu- 
larization), (c) NLO DGLAP evolution is repro- 
duced exactly at the inclusive level (c) MC per- 
forms NLO evolution by itself, using new exclu- 
sive NLO evolution kernels, without help of ex- 
ternal pretabulated inclusive PDFs. 

We refer the reader to Refs. [13114115116] for 
discussions of other aspects of the scheme of the 
MC modelling of the QCD NLO evolution not 
reported in this contribution. We would like also 
to point out to some similarities of this project to 
works in Refs. [17118] . 

2. NLO PSMC for C F part of NLO nons- 
inglet kernel 

We start from the ISR ladder of deep inelastic 
scattering of the lepton-hadron and our aim is to 
add 1st order corrections to the LO vertex in the 
middle of the ladder, see the Figure below 
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Figure 1. Deep inelastic lepton-proton scattering. 
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Figure 2. Correcting the last vertex in the LO 
ladder to NLO level. 



Let us start with correcting up to NLO level 
just one vertex (kernel) at the end of the ladder 
(next to hard process) in the LO PSMC, see the 
graph in Fig. [1] The multigluon distribution of 
the LO PSMC representing the inside of the LO 
PDF is 
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where Sisr is the Sudakov double log formfactor, 
the lightcone variable of the emitted gluon is at = 

k + k T 

2J-, the angular scale variable is <ij = -£■ and 
the actual distribution of the q — > Gq splitting is 

pf\ki) = 2C ; a5 fcT2^-, where l-z t =on/(l- 

Si=i a j)j an d Cf is colour factor. 

In Fig. [2] we illustrate again the situation in 
which the last (top) vertex is upgraded to NLO 
level. Hence in the leftmost graph this vertex 
includes already multiplicative virtual+soft NLO 
correction: 
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The two-gluon NLO correction coming from two 
graphs T , + T':< is defined as: 
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Here, the projection operator P, represented by 
the rectangular box in the graph i=i=i , , should 

be understood as operating on the amplitude 
squared, spin-summed: 
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This P operator is present in Fig.[2]in each ladder 
many times and its role is to simplify the exact 
matrix element to the LO approximation level all 
over the phase space. The role of the NLO cor- 
rection u[ is to undo the simplification done by 

P, just for 2 gluons. One could think that it is 
enough to include two leftmost graphs of Fig. [2j 
It is, however, necessary to include the sum over 
all LO spectators down to the beginning of the 
ladder, see Fig. [21 for the second gluon entering 

the NLO correction it' 

Altogether, the above NLO-corrected 
multigluon distribution reads as follows 
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where the MC weight component correcting the 
2-gluon distribution is: 



W(k 2) h) = 



1 9 i 1 9 

T-H \T"'"\ 



- 1 



and the virtual correction is /3n 
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Amazingly, in the expression for the PDF, 
D ( x ) = J27=al dLi P s nPn(ki,k 2 ,-kn)S(x - 
\Xj—i z j)> the integration of the NLO part 
J2j W(k n , kj) can be done analytically leading to: 
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recovering for the last emission precisely the NLO 
part (including virtuals) of standard DGLAP ker- 
nel y q q(u) defined according to: 
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In this way at the inclusive level the NLO stan- 
dard inclusive kernel of DGLAP is truly repro- 
duced for the last (n-th) vertex in the ladder. 

One can repeat the same procedure for any ver- 
tex number p\ in the middle of the ladder, sum 
up over pi , then apply the same for any two ver- 
tices p\ and pi and so on, such that finally all 
vertices in the ladder are at the NLO level. The 
above general procedure will yield the following 
multigluon distribution: 
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The above formula was already tested numeri- 
cally with three digit precision. 
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demonstrates 3-digit agreement, in units of the 
LO. 

Let us note that the Monte Carlo weight imple- 
menting NLO corrections is positive, with small 
dispersion and without any nasty tails. 

3. Adding CfCa part of the NLO nonsin- 
glet kernel of NLO PSMC 

Straightforward inclusion of the gluon pair dia- 
gram in the method of the previous section would 
be mathematically correct, but the Monte Carlo 
weight would be spoilt due to presence of Sudakov 
double logarithmic contribution +Sfsr in the 2- 
real-gluon correction 
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compensated by —Sfsr. component in the virtual 
correction 
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Resummation/exponentiation of FSR is therefore 
mandatory. Let us start again with the problem 
of upgrading to the NLO level the last vertex of 
the LO ladder. In the MC weight we add an extra 
term, which in the graphical form looks as follows 
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In the upper part of the figure above we 
present numerical results for the (inclusive) PDF, 
D(x,Q), from two Monte Carlo models, one 
with traditional inclusive evolution kernels and 
another one implementing our new NLO distri- 
bution of Eq. (JlJ . The trivial LO contribution is 
excluded from the comparison and the two curves 
represent the term in Eq. ([1]) with single W and 
with double W. Since the results cannot be dis- 
tinguished in the upper plot, we also show their 
ratio in the lower plot in the Figure above. The 
evolution runs from Qo = lOGeV to Q = ITeV, 
starting from D(x,Qo) = 5(1 — x). The ratio 



where the Sudakov term S FSR is subtracted in the 
virtual part: 
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f- = (i + 2K(A JSfl + v FSR -s FSR )) y* , 

and the FSR soft counterterm is subtracted in the 
2-real-gluon part: 
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It is remarkable that both above NLO correc- 
tions are free of any soft and collinear divergences, 



in spite of their complicated structure. This is 
thanks to a clever choice of the FSR counterterm 

!■{]■< , which is iterated in the FSR LO MC for 

the gluon emitted from the ladder. It was help- 
ful in the above exercise that we were employing 
angular ordering in the LO MC both for ISR and 
FSR parts. 

Altogether, the complete MC distribution, with 
the last LO vertex in the ladder upgraded to the 
NLO level looks as follows: 
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where NLO correction building blocks are 
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Again it is possible to check analytically, that the 
above apparently complicated MC weight repro- 
duces exactly the traditional NLO integrated ker- 
nels for the last vertex in the ladder. We have also 
performed numerical tests in which the above dis- 
tribution is modelled for the n = 1, 2 ISR gluons 
and any number of FSR gluons. The compari- 
son between MC and analytical result is in the 
following plot: 




It shows again the excellent agreement of the ex- 
clusive implementation of the NLO corrections 
in the PSMC and the analytical crosscheck, this 
time for the complete gluonstrahlung, including 
both Cp and CfCa parts. Again we have checked 
that the MC weight is positive and the distribu- 
tion of the MC weights is narrow. 

4. Summary and prospects 

Summarizing we state that the hrst serious fea- 
sibility study of the true NLO exclusive MC par- 
ton shower is almost complete for the non-singlet 
NLO DGLAP. The building block for the corre- 
sponding parton shower MC is tested numerically. 
Further work will cover the following areas: 

• Short range aim: Complete non-singlet, 
also for hard process. 

• Middle range aim: Complete singlet (quark- 
gluon transitions). 

• Optimize MC weight evaluation (CPU 
time). 



• Complete NLO MC for DIS at HERA and 
W/Z production at LHC. 

• Interface to standard NLO WS PDFs, 

• Extensions towards CCFM/BFKL. 
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